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I.  SUMMARY  OF  FINDINGS 

The  eftbrts  of  the  project  can  be  divided  into  three  separate  areas  of  endeavor. 
These  are  briefly  de.scribed  below.  It  is  intended  here  to  provide  enough  detail  to  give  the 
reader  a  grasp  of  the  main  thrust  of  the  work  without  getting  into  unnecessary  details. 

I..A.  DEVELOPMENT  OF  CONSTITUTIVE  THEORY 

The  behavior  of  a  material  must  be  consistent  with  restrictions  imposed  by  the 
second  law  and  the  principle  of  material  frame  difference.  If  is  the  Helmholtz  free 
energy,  the  entropy,  6  the  absolute  temperature,  T  the  second  Piola-Kirchoff  stress  tensor, 
and  E  the  Langrangian  strain  tensor,  the  second  law  may  be  written  as: 

Po  ^  -  Po  ^  ®')  ■  (l/9)F^q»  d  ^O  (1) 

F  is  the  deformation  gradient  and  is  the  initial  snow  density.  This  law  places  restrictions 
on  the  directions  in  which  processes  can  evolve.  Later  we  will  use  it  in  the  development 
of  the  constitutive  law.  The  principle  of  material  frame  indifference  is  the  statement  that 
the  properties  of  a  material  are  intrinsic,  i.e.  the  manner  in  which  it  responds  to  a  process 
(loading  or  deformation)  depends  on  the  material  itself  (Billington  and  Tate,  1981).  If  Q 
is  a  rotation  of  coordinates  from  one  coordinate  system  x  to  another  x*.  i.e. 

X*  =  Qx  (2) 

then  the  constitutive  law,  if  it  has  the  representation 

T  =  F  (E),  (3) 

must  transform  as 

T*  =  F(E*)  =  QF(E)Q^  (4) 

where  E*  is  the  strain  defined  relative  to  the  x*  coordinate  system.  Eqs.  1  and  4  need  tii 
be  satisfied  by  the  constitutive  law. 


We  begin  our  I'ormulation  by  assuming  ^  to  be  dependent  upon  the  stress  ’I',  the 
temperature  0,  and  an  internal  state  vector  This  state  vector, 

^  =  (?;.  ^2  •  ■  ■  ^n) 

is  a  set  of  scalar  valued  variables  which  describe  the  internal  state  of  the  material.  Later 
we  chose  the.se  to  repre.sent  variables  such  as  bond  diameter,  neck  length,  intergranular 
slip  distance,  etc.  Now  make  a  change  of  dependent  variables  by  introducing  the 
complimentary  energy  i/i,  which  has  the  value  T:E  -  4>  =  ip. 

Differentiating  vi  respect  to  time, 

V)*  -  (a^/aTT);!  +  (a0/ae)e  +  id4>/do  •  ^  (^’) 

and  substituting  into  the  .second  law  of  thermodynamics  allows  us  to  arrive  at  the  following 
results  if  the  second  law  is  to  be  satisfied: 


E  =  Po  ^V>/'3T  (7) 

T}  =-P(j  avi/ae  (8) 

Using  Equation  (7),  a  differential  change  in  the  strain  can  be  written  as: 
dE  =  Pg  (aVaTaT):  dT  +  p^  (aVacaT)dc  =  d^  E  +  d^E  (9) 

d^E  is  the  elastic  change  in  the  strain,  while  d^E  is  the  plastic  change  in  the  strain.  The 
compliance  tensor  M  is  defined  to  be 

M  =  p„(aVaTaT)  =  m  (c,e)  (lO) 

and  is  a  fourth  order  tensor  consisting  of  various  elastic  coefficients.  The  plastic  change 
can  be  broken  into  two  parts 

d^E  =  dM:T  +  dE^  (11) 

where  the  first  term  represents  a  strain  increment  due  to  changes  in  the  compliance  tensor, 
and  dE^  is  the  change  in  the  plastic  strain  (permanent  set).  If  the  material  compliance 
docs  not  change  significantly,  then  d^E  and  dE^  are  equal.  E^  is  the  strain  which  remains 
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«Iic’n  the  strcs.s  is  relc;isecl,  i.e.,  is  just  9,  ^). 

Now  define  the  vector  of  to  be  the  thermodynamic  conjugate  to  the  state  vector 
f  =  (12) 

Then,  one  can  show  that  (Hansen  and  Brown,  1988)  the  following  Maxwell 
relation  holds : 

attdT  -  aEldi  (13) 

Substituting  this  into  Eq.  9  and  integrating  gives,  after  some  algebra 
Po^  ~  E^:T  +  Pq0^  (14) 

where  is  a  stress  independent  term.  The  first  term  on  the  right  hand  side  is  the  elastic 
strain  energy  and  is  the  usual  quadratic  form  for  this  type  of  strain  energy.  The  form  of 
E  needs  to  he  determined.  This  is  best  approached  by  finding  a  specific  form  for  f.  One 
can  readily  determine  that 

f  =  p^  3xt>Jd^  +  1/2  T:(dM/d^):T  +  T:dE^/d^  (15) 

and  are  both  stress  independent  terms.  Therefore  the  above  equation  shows  that 
f  has  a  quadratic  dependence  on  T  if  M  is  a  function  of  as  it  surely  must  be,  and  that 
f  must  also  be  linear  in  T  as  required  by  the  last  term.  Given  Eq.  12,  ^  must  have  a  similar 
dependence  on  T. 

Using  the  above  relations,  and  requiring  to  be  a  scalar  valued  invariant  function 
of  T,  an  appropriate  form  for  is 

yj)  ^  +  (1/2)T:M:T  +  tr  (H  T)  (16) 

where  tr  (  )  is  the  trace  of  the  tensor  inside  the  parentheses.  H  is  a  second  order 

symmetric  tensor  which  changes  as  the  microstructurc  changes  during  inelastic  deformation, 
i.c.  it  is  a  function  of  ^  and  e. 

The  strain,  from  Eq.  7  becomes 
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M:T  +  H 


(17) 


E  =  Po 

and  we  see  that  H  is  closely  related  to  the  plastic  strain. 

The  clastic  strain  energy  part  of  V*  can  be  reduced  to  the  following  form 
1/2  T;M:T  =  {\+u)l2E  tr(T-)  -  i//2E  (trT)"  (IS) 

if  the  material  is  isotropic,  i/  and  E  are  Poisson’s  ratio  and  Young’s  modulus  and  arc 
functions  of  Let  Cj  and  be  the  coefficients  of  tr(T)"  and  (trT)"  cxpro.sscd  directly 


in  terms  of  elastic  propertic.s. 

Cj  =  1+V/2E,  C,  =  -V/2E  (19) 

They  can  be  expanded  as  a  Taylor  series  in 

^7  ~  ^  ^^77/  ■  ^i)  ^13i 

^2  ~  ^20  ^[^277  (^7  ■  ^7o)  ^2J7  (^7  '  ^70)"^]  (21) 

where  are  the  initial  values  of  each  state  variable.  In  the  above  an  odd  dependence 

on  was  assumed.  A  similar  form  for  11  is  assumed: 

^ij  =  ^ijo  (^7:  *  ^ijkSk  ^^k  '  (22) 


The  coefficients  in  the  above  equations  can  be  determined  from  experiments,  although  many 
will  be  assumed  to  be  negligible  based  on  physical  reasoning  and  symmetry  arguments. 
Finally,  evolution  equations  which  govern  the  rate  of  change  of  the  state  variables  ^  must 
be  found.  These  are  of  the  form 

^  =  ^(T,  e.  0  (2a) 

In  order  to  determine  these,  the  details  of  the  microstructure  of  the  material  must  be 
considered. 

MICROSTRUCTURAL  DESCRIPTION 
The  vector  ^  -  (^j<  4n)  consists  of  a  set  of  n  scalar  valued  variables  which 

describe  the  microstructurc  of  the  material.  For  the  purpicsc  of  this  study,  these  v.uiables  arc 


A 


assumed  to  be 


1  =  neck  length  vector 
r  =  neck  (bond)  radius 
A  =  intergranular  slip  distance  vector 
n^  =  bonds/grain  (coordination  number) 

These  variables  should  be  considered  as  mean  values  averaged  over  small  volume  elements. 
The  material  is  represented  as  a  collection  of  grains  joined  by  necks  of  a  definite  length  and 
radius.  When  two  ice  grains  are  brought  into  contact,  a  bond  forms  due  to  cohesion  procc.sses. 
This  neck  normally  grows  due  to  sintering  effects  and  evolves  into  a  necked  region  connecting  the 
two  grains.  Consequently  all  snow  that  has  been  on  the  ground  for  a  finite  time  has  this  type  of 
structure. 


(d)  'Jlress  componentG  dnd  (d)  npck  rtiGplact^mpnl 
unit  vectors  components 


Figure  1.  Schematic  of  typical  ice  grain  and  neck. 


Hcwcvcr,  the  ability  to  measure  bond  radius  and  neck  length  in  snow  is  still  a  difficult  task  and 
has  to  date  not  been  adequately  solved,  even  with  modern  image  analysis  techniques.  As  the 
theory  of  quantitative  stereology  and  the  capabilities  of  image  analysis  systems  continue  to  improve, 
these  problems  will  be  resolved. 

The  fundamental  geometry  of  an  ice  grain  and  a  neck  as  modelled  in  this  paper  is 
[)resented  in  Figure  1.  We  assume  most  of  the  deformation  takes  place  in  the  necks  which  are 
considerably  more  fiexible  than  the  grains. 

The  ice  grain  has  a  radius  R.  while  a  typical  neck  is  located  at  the  spherical  coordinates 
or,  3  and  has  a  radius  r  and  length  1.  Contrary  to  the  impression  given  in  the  figure,  the  neck 
length  is  typically  much  smaller  than  the  grain  radius.  The  stresses  acting  on  the  neck  have  the 
local  coordinate  components  aj^,  <72j,  and  ajj.  Using  an  appropriate  constitutive  relation  for  ice 
(Szyszkowski  and  Glockner,  1987)  the  rate  of  deformation  tensor  d  can  be  found  for  the  necked 
region.  TTic  full  details  of  this  constitutive  law  for  ice  cannot  be  given  here  due  to  space 
limitations.  A  more  general  constitutive  law  (Brown,  1987)  is  currently  under  development,  but 
the  formulation  of  Szyszkowski  and  Glochner  (1987)  has  been  shown  to  represent  the  properties 
of  ice  for  small  and  intermediate  strains.  The  rate  of  deformation  is  given  by  the  relation 


d  ^  d^"  +  d^  +  (T 

(24) 

where  d^  =  (1/E^)[(1 +v)a  -  /ilr(a)ll 

(25) 

d^  =  (1/V2)[‘^/(1-W))''(1  + 

(2b) 

d^  =  (l/v;)(a")" 

(27) 

P  r 

d  ,  d  ,  d  are  respectively  instantaneous  elastic,  plastic,  and  dehiyed  recoverable  (viscoelastic)  j),irts 
of  d.  a  is  the  slre.ss  tensor  in  the  neck,  of  which  and  are  the  only  nonzert^ 

Cv)rnponents.  E  and  v  are  elastic  properties  of  ice,  w  is  a  crack  damage  factor,  and  a  is  a  stress 
variable  in  the  neck.  Consult  Szyszkowski  and  Glockner  (ldS7)  for  a  mcme  detailed  discussion. 
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The  crack  damage  factor  w  characterizes  the  internal  damage  resulting  from  microcracKing. 

d  can  then  be  used  to  determine  the  deformation  of  the  neck,  i.e.  change  in  length, 
shearing  deformations,  and  change  in  radius  r.  The  motion  of  one  end  of  the  neck  relative 
to  other  end  is  given  hy  Xj,  and  A-)  give  the  shearing  motions  and  5  is  the  change  in  length 

of  the  neck.  Sec  Figure  1(a)  for  an  illustration  of  these  three  quantities. 

•  • 

The  rate  of  change  of  the  neck  length,  S,  and  the  two  intergranular  slip  rates,  Xj 

and  A-),  arc  related  to  d  in  the  following  manner: 

A,  =  2  d,jl,  i  =  1,2  (2S) 

In  arriving  at  these  relations,  it  was  assumed  that  rotational  effects  in  the  neck 
during  deformadon  is  negligible.  Normally  this  is  a  good  approximation  if  the  strains  in  the  neck 
arc  not  large.  It  should  be  mentioned  here  that  Xj  and  X2  represent  the  displacement  of  one  end 
of  a  neck  relative  to  the  other  end  in  a  direction  perpendicular  to  the  neck  axis.  Here  we  consider 
displacements  due  to  viscoplastic  deformation.  Slip  due  to  neck  fracture  and  glide  arc  not 
considered  here.  Studies  are  currently  underway  to  include  this  in  the  formulation. 

Thc.se  have  the  vector  forms  for  their  rates  of  change 

A  —  Xj  Hj  +  A-^a^  (29) 

5  =  6aj 

the  vectors  a^-  arc  local  orthogonal  unit  vectors  illustrated  in  Figure  1.  A  and  S  can  be 
uniquely  determined  from  the  tensor  d,  the  rate  of  deformation  tensor  for  the  neck. 

At  this  point  it  is  necessary  to  relate  the  neck  stress  a  to  the  globally  applied  stress 
T,  which  is  the  snow  stress.  Once  this  done,  the  formulation  is  essentially  complete.  To  this 
end,  we  define  a  probability  function  P(q,  n^)  which  gives  the  probability  that  the  point  on  the 
surlace  ol  a  grain  lies  on  a  point  where  a  neck  is  attached  to  the  grain.  I'his  function  must  be 
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(30) 


properly  normalized,  i.e. 

J  I  P(q,  /3.  nj)sin;3dQd/9  =  n(r/R)“ 

P  is  obviously  linear  in  nj.  For  the  purposes  of  this  study  we  will  assume  an  isotropic 
distribution  of  necks  around  the  grains,  so  that  P  becomes  a  constant  with  respect  to  q  and 
/3.  However,  it  is  worth  noting  that  anistropy  can  readily  be  modeled  by  making  P  variable 
in  a  and  p.  Indeed,  as  deformations  become  large,  P  would  normally  acquire  some 
anistropic  properties.  In  the  case  where  P  is  constant,  it  must  have  the  value  n^(r/2R)". 
If  T  is  the  stress  tensor  applied  to  the  snow,  the  stress  vector  applied  to  the  neck  can 
be  shown  to  be  (Hansen  and  Brown,  1986): 

=  a'  a,-  =  [a/p(a,^,nj)]n  T  (31) 

=  A  n  T 

n  is  the  unit  vector  parallel  to  the  direction  of  the  .neck  axis.  A  represents  the 
augmentation  factor  which  provides  for  the  increase  in  the  neck  slrc.sses  over  what  stresses 
are  distributed  over  the  snow,  a  is  the  density  ratio,  p,-/py  of  the  ice  density  and  snow 
density.  3^,  a%  and  a^  shown  in  Figure  2  have  the  values 
ay  =  cosasin^Cy  +  sinacos^e2-sina;0ej 

a,  =  -sinaCy  +  cosae,  (32) 

aj  =  cososin^Cy  +  sinasin^e,  +  cos^cj 

where  the  vectors  Cy  are  the  unit  cartesian  ba.se  vectors  for  the  x,  coordinates. 

0 

Under  a  three  dimensional  state  of  stres.s,  T  applied  to  the  snow,  the  strain  rate  F 
for  the  snow  can  be  calculated  using  the  above  equations.  First  on  a  neck  situated  at  the 
position  (a,/9)  can  be  calculated  by  using  Eq.  31.  From  thi.s,  the  rate  of  deroriiiation  tensor  d  Ibr 

the  neck  can  be  found  using  the  constitutive  equation  for  ice.  Eqs.  24-25.  Then  the  rate  of  change 

»  »  »  , 

of  the  microstructural  variables  A  and  1(1  =  6)  can  be  Ibund  with  Eqs.  28-29.  This  gives  the 


nitcs  A  and  S  in  the  global  cartesian  coordinate  system,  x-,  and  then  integrating  over  the  grain 

•  * 

surface  results  with  averaged  values  of  the  projections  of  A  and  1  in  the  global  coordinate 

directions.  Integration  over  an  octant  a<,  11/2,  Q  ^  p  <Y[J2  can  be  done  if  one  assumes  sufficient 

symmetry  of  the  probability  function  P(a.b).  These  averaged  values,  S  and  A  are 

m  n/2 

5  -  8((2P/p)V‘^3j  I  P(a,^)«(a,y3)sin;3dad^  (33) 

^  '^0 

A  =  ((2R/r)“/nj)  1  1  P(a,^)  A(a,/3)sin^dQd;9  (34) 

6  d 

•  * 

where  A  and  S  are  transformed  to  the  global  coordinate  components.  Once  this  is  done 
the  evolution  equations,  Eqs.  25,  are  essentially  determined.  With  these  determined,  the 
rates  of  change  of  the  compliance  tensor  (here  characterized  by  Cj  and  C,)  and  the  tensor 
H  can  be  determined  from  Eqs.  19-21,  and  then  Eq.  9  can  be  used  to  calculate  the  strain 
rate  E. 

The  above  full  program  has  not  yet  been  completed,  but  results  have  been  obtained 
for  simple  uniaxial  states  of  stress.  These  are  illustrated  in  the  following  .section. 

EVALUATION  FOR  UNIAXIAL  STRESS  STATES 
We  consider  here  the  strain  response  to  a  uniaxially  applied  load.  We  consider 
medium  density  snow  [p^  =  355  kg/m'^).  For  this  snow,  surface  section  analysis  has  determined 
the  following  initial  microstructural  variables 
R  =  0.484  mm 
r  =  0.23  mm 
1  =  0.19  mm 
nj  =  2.53 
a  =  2.58 

These  values  have  been  taken  from  a  study  (Man.sen  and  Brown,  1986)  utili/ing  an  image 
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analysis  system  to  evaluate  the  initial  microstructural  geometry  of  the  necks.  The  material 
coclTicients  for  the  constitutive  equation  for  the  icc  (Eqs.  24-27)  arc  provided  by  Szyszkowski  and 
Glockncr  (1986)  and  will  not  be  repeated  here.  Applying  the  previous  equations  to  the  case  of 
a  uniaxial  stress  state 

dEy^  =  2p^{C-)  +  C-))dTyy  -t-  2pp(3Cy/3^  -b  3C2/3^)  •  d^  +  dHyy  (^5) 

provides  for  an  increment  in  the  axial  strain  during  a  time  increment  when  the  stress  T j  y 
is  varied  or  the  microstructural  variables  are  varied. 

For  the  case  of  uniaxial  loading,  we  can  simplify  the  above  formulation  somewhat 
by  calculating  the  total  intergranular  slip  X  rather  than  the  two  components  Ay,  A-)  in  the  q  and 
p  directions  indicated  in  Figure  1.  Likewise,  rather  than  dealing  with  the  three  vector  components 
of  the  neck  length  vector  1,  we  can  deal  strictly  with  the  scalar  value,  1,  of  the  neck  length. 
Under  more  general  loading  conditions  where  deformations  in  all  their  coordinates  arc  of  concern, 
this  is  not  a  valid  simplification.  The  terms  aCy/SAandsC^/S^  then  require  the  forms 


dCjlaX  -  -  Cjjj  -b  3Cyjy(l  -  1^)^  (36) 

dC->/dX  ~  C-jyA  +  3C2jA( A-A^)“,  3C>/3l  =  ^3]]  'i' 

aH;y/3A  =  Hyy,  +  3Hy2^(A-A^)-  (38) 

3H;y/al  =  Hy;  -b  3Hy2y(l  -  IJ"  (29) 


Utilizing  experimentally  determined  values  of  the  coefficients  in  the  above  equations,  the 
rates  of  change  of  Hyy,  Cjj,  and  C2  can  be  obtained  as  the  microstructurc  changes,  since 
Hyy,  Cy,  and  C7  can  be  determined  directly  in  terms  of  Eqs.  36-39  and  the  rates  of  change 
of  X  and  1.  A  and  1  can  be  found  from  Eq.  28.  This  then  is  all  applied  to  Eq.  3.‘i  and 
integrated  numerically  with  respect  to  lime. 

Figure  2  illustrates  creep  curves  for  medium  density  snow  subjected  to  a  tensile 
strc.ss  of  0.01  MPa  for  a  prc.scribed  period  of  time  and  then  unloaded.  As  can  be  seen  the 
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material  exhibits  all  of  the  usual  characteristics  of  snow.  Also  contained  in  that  figure  is 
the  creep  curve  for  an  actual  experiment  carried  out  on  medium  density  snow. 
Unfortunately  the  data  was  collected  in  1974,  and  the  relevant  microstructural  variables  were 
not  recorded.  Consequently  it  is  not  known  if  this  data  is  relevant  to  the  case  used  here. 
Figure  3  shows  the  creep  respon.se  of  the  material.  The  response  cannot  be  compared  with 
available  data  in  a  definite  manner,  since  the  microstructure  of  the  snow  tested  was  not 
recorded.  The  data  at  MSU  was  acquired  in  a  testing  program  carried  out  between  1974 
and  1976.  While  the  creep  curves  for  that  data  was  similar  to  those  given  in  Figure  3,  the 
computed  strains  were  only  about  50%  of  the  measured  ones.  This  was  true  for  both  the 
axial  and  lateral  strains.  It  is  quite  po.ssible  the  snow  which  was  tested  had  incurred 
temperature  gradient  effects  and  was  poorly  bonded. 

The  constitutive  equation  still  needs  more  development  and  refinement  before  it  can 
acquire  a  truly  acceptable  form.  The  properties  of  snow  are  determined  in  part  by  the 
properties  of  the  matrix  material  ice.  Consequently  accurate  modelling  of  the  ice  properties 
is  needed.  However  the  snow  properties  are  also  strongly  inlluenced  by  the  geometry  of 
the  material  microstructure.  Correctly  modelling  this  is  the  more  difficult  of  the  two 
processes.  Work  is  currently  in  progress  to  obtain  a  more  complete  representation  of  the 
deformation  procc.sses  in  terms  of  what  happens  at  the  granular  level. 

This  type  of  constitutive  equation,  while  quite  involved,  can  be  used  to  more 
accurately  analyze  problems  involving  large  multitixial  deformations  of  snow.  Such  examples 
include  vehicle  mobility  over  snow  covered  terrain,  building  footing  settlement  into  snow 
cover,  and  impact  probIem.s.  With  ever  improving  computational  cajiability,  the  use  of 
microstructurally  based  stress  strain  equations  is  now  becoming  a  reality. 
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I.B.  EXPERIMENTAL  TECHNIQUE  FOR  ANALYZING  SNOW  MICROSTRUCTURE 

As  a  cohesive  granular  material  such  as  snow  is  deformed  under  the  influence  of 
applied  loads,  its  microstructure  changes.  In  fact,  even  as  the  material  is  at  rest  in  an 
unloaded  state,  its  microstructure  slowly  changes  as  the  material  attempts  to  lower  its  energy 
state  by  consolidating  to  reduce  its  surface  energy.  At  this  point  we  arc  primarily  interested 
in  the  changes  resulting  from  applied  loads.  As  loads  are  applied,  the  necks  connecting  the 
grains  are  strained,  their  lengths  and  diameters  change.  In  some  cases  they  fracture,  thereby 
initiating  intergranular  glide.  As  the  grains  are  pressed  together  under  compressive  loads, 
new  bonds  are  formed,  so  that  the  coordination  number  (average  bonds/grain)  changes. 
Under  certain  conditions  grains  can  fracture,  thereby  causing  changes  in  mean  grain  size. 

The  constitutive  law  discussed  in  Section  LA  is  described  in  terms  of  these 
microstructural  variables.  As  a  consequence,  we  needed  to  develop  a  means  of 
experimentally  determining  the  variation  of  these  variables  as  a  result  of  load  induced 
deformations.  This  is  necessary  if  we  were  able  to  determine  if  the  evolution  equations 
(Eqs.  23)  were  accurately  predicting  these  microstructural  changes.  Tliis  section  describes 
in  some  detail  the  method  we  developed  to  achieve  these  ends. 

The  internal-state  variables  must  de.scribe  the  current  microstructural  state  of  the 
material  and  be  capable  of  representing  average  measures  o  the  structural  rearrangements 
taking  place  within  the  snow  cover.  The  variable  selection  procc.ss  attempted  to  include 
the  most  significant  aspects  of  the  granular  structure  and  the  dynamic  environment  of  the 
material  under  high-rate  deformation. 

The  variables  cho.sen  should  be  able  to  characterize  the  dominant  deformation 
mechanisms  such  as  bond  fracture,  intergranular  glide,  and  neck  growth  at  the  bonds. 
Furthermore,  they  must  account  for  various  phenomena  known  to  occur  under  compressive 


loadings  such  as: 

(1)  The  effects  of  pore  pressure, 

(2)  A  locking  mechanism  for  grains  during  intergranular  glide, 

(3)  Coupling  of  the  deviatoric  and  volumetric  responses, 

(4)  Work  hardening,  and 

(5)  Local  inertial  effects. 

Finally,  the  state  variables  must  be  chosen  so  as  to  make  measuring  them  a  feasible  task. 

Before  proceeding  with  the  selection  of  the  state  variables,  it  is  necessary  to  clarity 
the  notation  to  be  used.  When  discussing  the  behavior  of  dry  snow  in  the  general 
macroscopic  sense,  the  material  is  considered  to  be  a  multiphase  mixture  consisting  of  an 
ice  phase  and  an  air  phase.  Each  phase  is  denoted  by  a  subscript  as  follows:  a,  air  phase; 
i,  ice  phase.  This  convention  is  in  accordance  with  that  used  in  mixture-theory  studies  of 
mechanical  properties  of  snow. 

An  alternative  approach  will  be  used  when  measuring  the  internal-state  variables  in 
the  snow  at  the  granular  level.  Here,  the  snow  is  considered  to  be  a  three-phase  material 
consisting  of  an  ice-grain  phase,  ice-neck  phase,  and  air  phase.  The  three  phases  arc 
denoted  by  the  following  Greek  subscripts:  a,  ice -grain  phase;  p,  ice-neck  phase;  and  7, 
air  phase. 

In  the  following  discussions,  no  attempt  is  made  to  state  which  of  the  two  previous 
approaches  is  being  used  since  the  subscripts  arc  self-explanatory.  Finally,  in  quantitative 
strcology,  it  is  nccc.ssary  to  differentiate  between  length,  area,  and  volume  measurements. 
Therefore,  let  the  following  capitalized  subscripts  denote  this  difference:  L,  length 
measurement;  A,  area  measurement;  and  v,  volume  measurement. 
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Taking  into  account  the  desired  properties  of  the  state  variables,  the  following 
parameters  were  chosen  to  formulate  the  statistical  model  for  snow.  A  discussion  of  their 
significance  and  measurability  is  given  later. 

a  =  p^lp  —  density  ratio  where  p^  =  the  density  of  ice,  and  p-  the  density  of  snow, 
r  mean  bond  radius , 

h  mean  bond  length, 

7  mean  intergranular  slip  distance, 

L  mean  intercept  length, 
nj  mean  number  of  bonds  per  grain, 

mean  number  of  grains  per  unit  volume 
V  mean  volume  of  a  single  grain,  and 
S  mean  surface  area  per  unit  volume. 

MEASURING  THE  STATE  VARIABLES 

Before  discussing  each  internal-state  variable,  it  is  necessary  to 
describe  the  streological  measurements  which  are  required  to  obtain  the 
desired  information.  As  mentioned  previously,  snow  at  the  granular  level 
is  considered  to  be  a  three-phase  material  consisting  of  ice  grains,  ice 
necks,  and  air.  Kry  (1975)  developed  an  operational  definition  to  identify 
grain  bonds  in  a  two-dimensional  surface  section.  If  a  grain  bond  is  cut 
by  a  section  plane,  it  will  appear  as  a  line  connecting  opposite  edges  of 
the  ice.  Three  criteria  are  required  to  identify  these  bonds: 

(1)  A  minimum  constriction  must  exist;  a  30Z  constriction  on  the  plane 
is  used  as  a  cut-off. 

(2)  Both  edges  of  the  ice  must  show  the  constriction. 

(3)  The  notches  on  each  edge  must  point  approximately  towards  each  other. 
The  neck  region  of  a  grain  bond  is  defined  as  that  area  surrounding  the 
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Figure  4.  Typical  screen  image  analysis  of  snow  surface  section. 
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bond  where  the  grain  goes  from  convex  to  concave  with  respect  to  the 
outward  normal.  The  bond  is  defined  as  the  region  of  minimum  constriction. 
Figure  4  shows  an  image -enhanced  snow  sample  which  delineates  the  three 
phases  of  interest.  Using  the  enhanced  image,  the  following  measurements 
can  now  be  made.  The  code  contained  within  the  parentheses  indicates  the 
method  of  analysis,  e.g.  visual  count  (VC)  or  computerized  image  analysis 
(lA)  . 

Area  fraction  of  points  falling  in  the  ice-grain  phase  (lA). 

Area  fraction  of  points  falling  in  the  ice-neck  phase  (lA). 

P.,  Area  fraction  of  points  falling  in  the  air  phase  (lA). 

E  =  Harmonic  mean  of  lines  representing  grain  bonds  in  the  surface  section  (lA). 

Hence,  62  is  defined  as  the  two-dimensional  bond  diameter. 

Number  of  interceptions  of  ice  grains  in  the  microstructure  per  unit  length  of  a 
random  test  line  (lA)  (See  Fig.  5). 

Number  of  grains  in  the  test  area  (VC) 

N^  Number  of  bonds  in  the  test  area  (VC). 

f2(n)  Probability  distribution  of  the  number  of  bonds  per  grain  cut  by  the  section  plane 
(VC). 

It  should  be  noted  that  all  mca.suremcnts  which  require  a  visual  count  can  be  performed 
on  the  order  of  minutes.  The  time-consuming  part  of  the  analysis  will  be  in  the  area  of 
image  enhancement. 

At  this  point  one  can  proceed  with  a  discussion  of  each  of  the  stale  variables. 
Density  Ratio  (a) 

The  density  ratio  is  defined  by  the  following  expre.s.sion: 

a  =  p^lp  (40) 
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The  density  of  snow  has  been  the  primary  measure  of  the  "state"  of  the  material  in  many 
of  the  past  investigations  of  the  mechanical  behavior  of  snow.  This  is  primarily  because  of 
the  ease  and  quickness  of  density  measurement  in  the  field.  However,  recent  papers 
(Gubler,  1978)  have  shown  that  density  alone  cannot  adequately  describe  the  state  of  the 
material  under  certain  conditions.  This  is  especially  true  at  low  densities  where  the  impact 
of  other  state  variables  is  more  significant.  However,  the  importance  of  density  increa.scs 
as  a  approaches  unity. 

The  value  of  the  density  ratio  obtained  in  the  field  may  be  cross  checked  using  the 
following  streological  relationship; 

a  =  1/(P,  +  P,)  (41) 

This  equation  neglects  the  mass  of  the  air  phase. 

Mean  bond  radius  (r) 

The  mean  three-dimensional  bond  radius  is  a  significant  parameter  of  pressure 
sintering  as  well  as  bond  strength  and  resistance  to  fracture.  Brown  (1979)  used  this 
parameter  in  a  volumetric  constitutive  law  based  on  neck  growth.  Pullman  (1953)  derived 
the  necessary  relations  with  the  assumption  that  the  grain  bonds  are  circular  disks  (1975) 
has  shown  that  this  idealization  yields  self-consistent  results. 

The  governing  equation  for  the  mean  bond  radius  as  derived  by  Pullman  is 
r  =  II/(4E)  (42) 

Mcun  bond  length 

The  mean  bond  length  is  al.so  a  significant  parameter  for  modeling  pressure  sintering. 
This  was  used  in  Brown’s  "neck-growth"  volumetric  constitutive  law.  However,  until  now 
the  authors  are  unaware  of  any  attempts  to  evaluate  this  parameter  stereologically.  The 
statistical  model  of  bond  length  is  developed  by  idealizing  the  necks  of  grains  to  be 
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cylinders.  This  is  based  on  the  fact  that  self-consistent  results  have  been  obtained  by 
idealizing  the  grain  bonds  as  disks.  Clearly  the  neck  regions  of  an  ice  grain  arc  not 
perfectly  cylindrical.  However,  the  definition  still  provides  a  meaningful  relationship  when 
comparing  differing  types  of  snow  in  that  it  is  a  consistent  measure  of  bond  length. 
Pullman  (1953)  has  derived  a  relationship  for  the  number  of  bonds  per  unit  volume  ba.scd 
on  knowledge  of  the  number  of  bonds  per  unit  area.  This  relationship  is 
N,v  =  8EN,^/^“  (43) 

where  is  the  number  of  bonds  per  unit  area.  Underwood  (1970)  has  shown  that  a 
statistically  exact  equation  for  determining  the  volume  fraction  of  the  ath  constituent  is 
given  by 

=  P.  (44) 

Therefore,  by  introducing  the  neck  region  of  an  ice  grain  as  a  separate  pha.se,  one  can 
determine  the  volume  fraction  of  the  necks.  The  idealized  bond  length  can  then  be 
determined  by  the  following  expression: 

h  =  P,/(N,,nr)  (45) 

Mean  intergranular  slip  distance  (A) 

The  mean  intergranular  slip  distance  is  of  major  importance  for 
characterizing  intergranular  glide.  As  this  distance  decreases,  pressure 
sintering  becomes  more  dominant  in  the  deformation  process.  Furthermore, 
this  parameter  must  be  inherently  related  to  the  locking  phenomenon  of 
intergranular  glide.  The  stereological  relationship  as  defined  by  Fullman 
(1953)  is  given  by 

A  =  (a-  l)/(aN^)  (46) 

This  equation  is  valid  regardless  of  size,  shape,  or  the  respective  distributions  of  the  grains. 

Mean  intercept  length  (L) 
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The  mean  intercept  length  provides  a  measure  of  the  grain-size  and,  as  a  result, 
complements  the  mean  intergranular  slip  distance.  The  combination  of  these  two 
parameters  should  provide  a  physical  basis  for  determining  a  critical  "locking"  density  for 
grains  during  intergranular  glide.  Furthermore,  pressure  sintering  of  neck  regions  w’ill  cause 
changes  in  L. 

Tlie  mean  intercept  length  as  derived  by  Underwood  (1970)  is  given  by 

L  =  (P,,  +  P,)N^  (47) 

Again,  this  equation  makes  no  assumptions  about  the  shape,  size,  or  distribution  of  tr.c 
grains. 

Remaining  state  variables 

The  remaining  internal-state  variables  must  be  discu.ssed  together  because  of  their 
close  dependence  on  each  other.  The  variables  are  given  by 

(i)  Number  of  grains  per  unit  volume  (N^y) 

(ii)  Mean  number  of  bonds  per  grain  (n^) 

(iii)  Mean  volume  of  a  single  grain  (V) 

(iv)  Mean  surface  area  per  unit  volume  (S). 

The  mean  number  of  bonds  per  grain  is  a  strong  measure  of  the  degree  of  grain 
mobility  and  fracture  strength.  The  number  of  grains  per  unit  volume  and  the  mean  volume 
of  a  grain  are  important  parameters  for  modeling  intergrai.ular  glide  including  local  inertial 
effects.  Finally,  the  mean  sui  ace  area  per  unit  volume  is  extremely  important  for 
computing  energy  absorption  by  snow  during  high-rate  deformation.  Under  high  strain- 
rates,  the  snow  will  undergo  brittle  fracture  at  the  grain  bonds.  This  ctiuscs  a  signi'icant 
ii.erease  in  the  surface  area  which  in  turn  is  directly  related  to  the  free  energy. 
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Until  recently,  no  method  was  available  for  determining  the  mean  number  of  bonds 
per  grain  from  a  section  plane.  Gubler  (1978)  derived  a  technique  for  determining  this 
value  based  on  comparing  a  theoretical  distribution  of  the  number  of  bonds  per  grain,  r^(n), 
seen  in  a  section  plane  with  the  actual  distribution.  This  paper  appears  to  be  a  pioneering 
effort  in  this  area.  The  technique  involved  calculating  the  probability,  p,  that  if  a  grain  with 
a  coordination  number  n  =  1  is  cut  by  a  section  plane,  its  bond  also  appears  in  the  section 
plane.  The  bond  areas  are  assumed  to  be  random  and  isotropically  distributed  o;i  the  grain 
surface.  Furthermore,  the  bonds  are  considered  to  be  small  compared  to  the  grain  surface. 
This  is  consistent  with  Kry’s  (1975)  definition  of  a  bond. 

Unfortunately,  Gubler’s  approach  contains  one  significant  weakness  in  that  the 
stercological  relationship  used  to  obtain  the  number  of  grains  per  unit  volume  is  based  on 
all  grains  having  the  same  size  and  shape.  This  is  hardly  the  case  for  alpine  snow,  as  highly 
faceted  crystals  as  well  as  rounded  grains  appear  throughout  the  snow  cover  in  many  sizes. 

The  problem  of  accurately  determining  the  number  of  grains  per  unit  volume  from 
surface  sections  for  a  collection  of  particles  of  arbitrary  size  and  shape  is  extremely  difficult. 
Most  authors  have  a.ssumed  the  particles  to  be  all  one  size  and  shape  for  which  many 
solutions  arc  known.  Some  solutions  arc  also  known  for  particles  of  one  shape  which  obey 
a  log-normal  distribution  in  size.  DcHoff  (1964-1965)  has  provided  solutions  for  this  ca.se 
for  .some  simple  .shape.s.  Finally,  Hilliard  (1968)  has  developed  a  technique  for  determining 
N^^.  for  an  arbitrary  collection  of  sizes  and  shapes  provided  the  relative  frequency  of  the 
various  shapes  is  known.  Unfortunately,  this  theory  is  unable  to  produce  specific  results 
since  shape  factors  nece.ssary  to  apply  the  theory  have  not  been  developed. 

The  approach  taken  here  to  determine  will  avoid  any  specific  identification  of 
grain  shape  other  than  to  obtain  a  first  gue.ss,  as  this  is  virtually  impossible  for  alpine  snow. 
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The  formulation  follows  Gublcr’s  approach  with  some  notable  exceptions. 


To  begin,  the  probability  of  finding  a  grain  with  a  coordination  number  between  n 
and  n  +  dn  as  derived  by  Gubler  is  given  by 

(48) 

where  n  is  the  three-dimensional  coordination  number,  N,  is  a  normalizing  constant,  and 
i  is  a  free  parameter  of  the  distribution. 


fj(n)  =  Njn'^'^exp(-^n-') 


and 


^  - - 1 

n  n_j 


3  1 
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The  normalizing  constant  may  be  determined  numerically  by  requiring 
12 


(49) 


n  =  l 


This  is  based  on  the  premise  that  each  ice  grain  in  a  snow  cover  has  at  least  one  bond  and 
at  most  12  bonds,  as  in  the  case  for  ice. 

The  resulting  two-dimensional  probability  distribution  as  derived  by  Gubler  is 

12 

f,(l)  = 

k=l 

where  p  is  defined  as  the  probability  that,  if  a  grain  with  a  coordination  number  n 
=  1  is  cut  by  a  section,  its  bond  also  appears  in  the  section,  and  N->  is  a  normalizing 
constant.  The  normalizing  constant  N7  does  not  appear  in  Gublcr’s  work.  The 
rea.son  for  inserting  this  is  as  follows.  Equation  (50)  without  N7  represents  a 
binomial  distribution  of  the  probability  p  multiplied  by  the  three-dimensional 
probability  distribution  fj(k).  Therefore,  the  following  equality  holds: 

12 

H  hO)  =  J  (51) 

1  =  0 


UJ 


P(1  -  p)^-^f3(k) 


(50) 


00 


However,  when  cutting  a  section  of  ice  with  grains  having  a  coordination  number  of  12,  it 
is  physically  possible  to  cut  at  most  six  bonds.  Hence,  one  can  require 
6 

y  f2(l)  =  1  (52) 

The  probability  p  is  given  by  the  following  (Hilliard,  1968): 

P  =  sA  (53) 

where 

=  fl-N^b/N^ 

=  4n^/N^ 
and 

b  =  2(r^  -  (dy2)^)^^' 

Notice  that  p  is  a  function  of  the  number  of  grains  per  unit  volume,  It  is  at  this  point 
where  all  previous  investigations  assume  specific  shape  and  size  distributions  and  where  this 
theory  deviates. 

Underwood  (1970)  provided  a  table  for  determining  for  particles  of  one  size  and  shape. 
TTiis  is  the  approach  followed  by  Gublcr.  The  general  equation  is  given  by 

N<.v  =  C(N-’,^/N^)  (54) 

where  C  is  a  coefficient  depending  only  on  shape.  It  is  interesting  to  note  that  C  takes  on 
a  fairly  narrow  range  of  values  for  a  wide  variety  of  shapes,  typically  ranging  from  0.4  to 
0.8  for  shapes  pertinent  to  alpine  snow.  Therefore,  one  can  use  Equation  (54)  to  obtain 
a  first  estimate  of  N^j,.  This  in  turn  fixes  the  probability,  p,  as  defined  in  Gubler’s  paper. 

Now  recall  an  earlier  relationship  derived  by  Pullman  (1953)  for  the  number  of 
bonds  per  unit  volume 

N,,  =  8ENp^/n-  (55) 
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The  above  equation  is  based  on  the  assumption  that  the  bonds  may  be  considered  to  be 
a  polydispersed  system  of  thin  disks.  This  is  precisely  what  has  been  assumed  throughout 
this  paper.  Using  this  result,  the  mean  number  of  bonds  per  grain  is  then  Oxed  by  the 
equation 

n^  =  2N,,/N,v  (56) 

By  neglecting  this  equation  and  allowing  to  vary,  Gubler  hac  an  underconstrained  system. 
This  could  lead  to  uniqueness  problems  and  produce  erroneous  results. 

By  invoking  Equation  (55),  the  parameters  and  i  are  varied  rather  than  nj  and 
i  to  determine  the  proper  theoretical  density  distribution  f^Cb).  Hence,  a  method  has  been 
developed  for  determining  N^^,  without  specific  consideration  of  the  individual  grain  shapes 
and  sizes  other  than  to  obtain  a  first  guess. 

Once  is  determined,  the  mean  volume  of  a  single  grain  including  the  "necks"  is  given 
by 

V  =  (P,  +  P,)/N,v  (57) 

The  neck  regions  in  this  variable  are  included  since  the  neck  is  assumed  to  become  an 
integral  part  of  the  grain  when  a  bond  fractures.  Furthermore,  the  neck  volume  is  typically 
on  the  order  of  5-10%  of  the  total  grain  volume.  Therefore,  any  errors  caused  by  this 
assumption  are  small. 

The  mean  surface  area  per  unit  volume  can  be  derived  by  considering  each  of  the 
grains  as  being  detached  from  its  neighbors  and  then  subtracting  from  the  surface  area 
contained  by  the  grain  bond.s.  Underwood  (1970)  derived  an  exprc.ssion  for  the  surface  area 
for  a  .system  of  detached  grains  given  by 

S  =  4N^  (58) 
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Therefore,  subtracting  the  area  contained  by  the  grain  bonds  from  the  above  expression 
gives  the  correct  surface  area  per  unit  volume 

S  =  4N<^  -  2n-b3N„v  (59) 

application 

The  theory  presented  above  has  been  written  into  a  computer  code  which  uses  a 
high-speed  image  analyzer  located  at  Montana  State  University.  This  image  analyzer  is 
capable  of  digitizing  a  photograph  of  a  section  into  approximately  300,000  points  and 
assigning  each  point  one  of  256  possible  gray  levels.  The  digitized  section  is  then  image- 
enhanced  and  stored  on  a  floppy  disk  for  future  analysis. 

By  using  the  image  analyzer,  exact  point  counts  and  extremely  accurate  length 
measurements  can  be  made  very  quickly.  This  allows  large  areas  of  a  section  containing 
numerous  grains  to  be  analyzed  which  increases  the  statistical  accuracy. 

The  following  data  are  an  example  of  the  results  of  the  theory  using  the  image 
analyzer.  A  surface  section  was  taken  from  a  box  of  alpine  snow  stored  in  a  cold  room  at  - 
12^C  over  the  summer  at  Montana  State  University.  Two  areas  of  the  section  were  then 
analyzed  using  the  image  analyzer.  The  areas  of  analysis  are  shown  by  the  dashed  and  solid 
lines  in  Figure  4,  respectively.  The  two  areas  overlap  and  hence  are  not  statistically 
independent.  However,  the  large  area  does  provide  a  check  on  the  statistical  accuracy  of 
the  smaller  area.  The  data  taken  from  the  image  analysis  of  the  two  areas  are  shown 
below. 


Area  1 

Area  2 

Area  of  analysis 

64.4  mm“ 

133.9  mm" 

Ice-grain  area  fraction 

0.350 

0.363 

Ice-neck  fraction 

0.038 

0.043 
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Air-phase  area  fraction  0.612  0.594 

Mean  two-dimensional  bond  length  0.307  mm  0.331  mm 

Harmonic  mean  two-dimensional 

bond  length  3.42  mm'^  3.30  mm'^ 

Number  of  grains  per  unit  length 

of  a  random  test  line  0.803  mm'-^  0.776  mm'^ 

The  data  generated  by  the  image  analysis  of  the  surface  section  are  in  good 
agreement  for  both  areas  analyzed.  This  indicates  that  the  smaller  area  was  representative 
of  the  overall  properties  of  the  section. 

Using  the  above  information,  the  two-dimensional  theoretical  probability  distribution 
for  grain  coordination  number  can  be  computed.  The  results  for  each  area  analyzed  along 
with  the  measured  distributions  are  shown  below. 


Area  1 

No.  of 

Measured 

Theoretical 

Bonds 

distribution 

distribution 

n  =  0 

0.518 

0.515 

n  =  1 

0.291 

0.313 

n  =  2 

0.139 

0.102 

n  >  3 

0.052 

0.070 

Area  2 

No.  of 

Measured 

Theoretical 

Bonds 

distribution 

distribution 

n  =  0 

0.490 

0.489 
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n  =  1 

0.327 

0.327 

n  =  2 

0.112 

0.111 

n  >  3 

0.071 

0.073 

The  theoretical  distributions  for  the  number  of  bonds  per  grain,  as  seen  in  the 
section  plane,  compared  very  favorably  with  the  measured  distributions.  For  the  larger 
area  of  analysis,  the  error  is  approximately  1%  for  grains  with  0,  1,  or  2  bonds,  respectively. 
These  grains  account  for  approximately  93%  of  the  grains  in  the  section  analyzed. 

Using  the  theoretical  data  generated  above  in  conjunction  with  the  previous 
information  determined  by  image  analysis,  the  internal-state  variables  may  be  readily 
computed.  These  are  given  by  the  following: 


Area  1 

Area  2 

Density  ratio 

2.58 

2.46 

Mean  three-dimensional  bond  radius 

0.230  mm 

0.238  mm 

Mean  three-dimensional  bond  length 

0.193  mm 

0.182  mm 

Mean  intergranular  slip  distance 

0.762  mm 

0.766  mm 

Mean  intercept  length 

0.484  mm 

0.523  mm 

Mean  number  of  bonds  per  grain 

2.04 

2.26 

Mean  number  of  grains  per 

unit  volume 

1.16  mm'"^ 

1.17  mm'-^ 

Mean  grain  volume 

0.334  mm-^ 

0.345  mm^ 

Mean  surface  area  per  unit  volume 

2.43  mm'^ 

2.16  mm'^ 

The  repeatability  of  the  data  for  the  internal-state  variables  for  the  two  areas 
analyzed  is  quite  strong  with  the  difference  between  the  two  typically  being  on  the  order 
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of  10%  or  less.  Also,  the  density  ratio  determined  by  a  field  measurement  of  this  particular 
snow  sample  was  2.29.  This  differs  from  the  stereological  measurements  by  13%  ar.d  7%, 
respectively  for  areas  1  and  2.  It  is  not  surprising  that  the  field-measured  density  was 
slightly  higher  than  that  determined  by  image  analysis,  since  some  compaction  of  the  snow 
will  occur  when  taking  field  data.  Furthermore,  etching  of  the  ice  grains  during  sectioning 
will  cause  the  image  analysis  to  predict  a  slightly  lower  density  than  the  actual  value. 

The  above  technique  was  then  used  to  measure  changes  in  the  microstructural 
parameters  during  deformation.  It  was  felt  that  the  most  substantial  changes  in 
microstructure  could  be  achieved  by  means  of  confined  compression  tests.  A  scries  of  tests 
were  run  to  compress  snow  from  an  initial  density  of  about  330  kg/m-^  to  as  high  as  650 
kg/m-^.  TTiis  upper  value  represented  the  densification  achieved  when  the  load  capability 
(10,000  lb)  of  the  testing  machine  was  reached.  In  one  test  the  deformation  was  stopped 
at  2000  pounds.  In  other  tests  the  loading  was  stopped  at  5000  lb  and  at  10,000  lb.  In 
order  to  obtain  statistically  viable  data,  five  tests  were  made  for  each  m:iximum  load. 

Prior  to  each  test,  a  surface  section  was  made  to  determine  the  initial 
microstructural  parameters.  After  each  test,  the  sample  was  cut  in  half  and  samplc.s 
removed  for  surface  sections.  These  surface  sections  were  then  image  analyzed  to 
determine  the  variation  in  the  microstructure  due  to  metamorphism.  For  example.  Figures 
5  and  6  demonstrate  the  variation  of  the  coordination  number  and  neck  length  with  density. 
TTic.se  are  typical  of  results  achieved.  A  more  detailed  descriptive  of  all  microstructural 
variables  is  forthcoming  in  future  publications. 

III.  C  Heat  and  Mass  Transport  in  Snow 

As  was  indicated  earlier,  the  microstructurc  of  snow  undergoes  changes  even  when 
it  is  unstrc,ssed.  Again,  since  microstructurc  determines  the  mechanical  properties,  the 
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DENSITY  (g/cc) 


Figure  5.  Variation  of  coordination  number  n3  with  density  during  compression 
of  snow. 


question  of  which  heat  and  mass  transfer  processes  effect  the  properties  was  of  concern. 
Consequently,  it  was  decided  to  adapt  a  mixture  theoi'y  to  investigate  this.  The  detaik  of  tlx; 
theoretical  basis  behind  the  theory  are  too  involved  to  present  here.  The  theory  is  very 
involved  and  a  considerable  development  effort  was  required.  Some  of  this  was  done  under 
sponsorship  of  an  NSF  grant,  and  the  application  of  this  theory  to  snow  was  done  under 
this  ARO  grant.  For  a  more  detailed  development  of  the  theory,  the  reader  is  referred  to 
the  papers  by  Adams  and  Brown  (1988  a,b). 

The  material  is  assumed  to  consist  of  a  mixture  of  constituents,  in  this  case  the 
vapor  phase  and  the  solid  phase,  ice.  For  each  constituent,  the  balance  equations  for  mass, 
momentum,  and  energy  are  postulated.  A  constitutive  behavior  for  each  constituent  is 
assumed,  and  the  second  law  of  thermodynamics  is  then  used  to  place  restrictions  on  these 
constitutive  relations.  The  conditions  of  thermomechanical  equilibrium  and  stable 
equilibrium  are  also  investigated  to  further  restrict  these  relations.  The  resulting  set  of 
constitutive  equations  and  balance  principles  form  a  set  of  coupled  nonlinear  differential 
equations  which  must  be  solved  numerically. 

This  theory  was  then  applied  to  three  sets  of  problems: 

(1)  temperature  gradient  effects 

(2)  density  gradient  effects 

(3)  fabric  layering  effects 

The  first  one,  often  referred  to  as  temperature  gradient  metamorphism.  A  temperature 
gradient  cau.ses  a  flux  of  mass  and  heat  in  the  direction  of  the  negative  temperature 
gradient.  This  produces  a  rectystallization  of  the  snow,  of  ten  resulting  with  substantial 
reductions  in  strength.  This  would  have  a  negative  effect  on  mobility  of  military  vehicles 
over  snow  covered  terrain. 
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The  other  two  processes,  which  we  refer  to  as  density  gradient  mctamorphisni  and 
fabric  gradient  metamorphism  result  due  to  variations  of  material  properties 
(inhomogeneities)  in  the  snow  cover.  Density  gradients  in  the  snow  cover  causes  mass  to 
move  from  the  lower  density  snow  to  the  higher  density  snow  as  the  material  attempts  to 
lower  its  surface  energy.  In  the  last  process,  vapor  is  transported  from  fine  grained  snow' 
to  course  grained  snow,  since  fine  grained  snow  has  a  smaller  ■  adius  of  curvature  and  hence 
a  higher  surface  energy. 

All  three  of  these  processes  as  well  as  others  contribute  simultaneously  to  the 
overall  flux  of  heat  and  mass  through  the  snow,  with  subsequent  alterations  of  the 
mechanical  properties.  In  an  attempt  to  learn  more  about  these  processes  and  to  better 
quantify  these,  the  mixture  theory  discussed  above  was  specialized  to  the  case  of  snow  and 
used  to  evaluate  each  process.  These  arc  briefly  discu.ssed  below. 

Temperature  Gradient  Metamorphi-sm: 

In  this  case,  a  one  meter  deep  snow  cover  of  homogeneous  snow  was  suddenly  cold 
from  1273^C  to  25.3^C  on  the  top  surface  while  the  bottom  was  maintained  at  27.V’C. 
Figure  7  shows  the  time  dependent  response  of  the  temperature  as  a  function  of  lime, 
while  Figures  8  and  9  show  respectively  the  time  dependent  vapor  fiux  and  condensation 
rate  (rate  of  increase  of  volume  of  ice)  in  the  snow  cover.  The  pictures  show  vividly  how 
the  establishment  of  a  temperature  gradient  caascs  a  flux  and  subsequent  rccondcnsation 
of  vapor  throughout  the  snow  cover.  Although  not  shown  here,  the  theory  also  brings  out 
other  more  subtle  results.  For  instance,  when  gravity  is  included  in  the  calculations,  the 
effect  of  the  stre.sscs  cau.sed  by  overburden  loads  increases  the  surface  encigy  slightly  in  the 
snow  near  the  bottom.  This  itself  causes  sublimation  of  vapor  off  the  crystals  near  the 
bottom  and  a  migration  of  vapor  upward. 
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Figure  9.  Time  dependent  response  of  condensation  rate  in  layered  snow  cover 
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Figure  10.  Density  profile  in  layered  snow  cover 
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Density  Gradient  Metamorphism: 

This  effect  has  not  been  recognized  until  this  study  was  implemented.  Less  dense 
snow  has  more  free  surface  area  in  it  than  dense  snow,  thereby  increasing  the  surface 
energy  for  the  Ic.ss  dense  snow.  In  order  for  a  material  to  decrease  its  free  energy,  it 
naturally  tries  to  decrease  its  surface  energy.  In  the  case  of  snow,  this  causes  a  flu.\  of 
vapor  from  low  density  regions  to  high  density  regions.  Given  enough  time,  snow  would 
eventually  consolidate  to  ice. 

In  order  to  study  this  effect,  the  layered  snow  cover  shown  in  Figure  10  was 
analyzed  using  the  mixture  theory.  The  results  shown  in  Figure  11  shows  that  over  time 
there  is  a  definite  migration  of  mass  from  the  less  dense  layers  into  the  dense  snow.  This 
causes  a  mass  loss  in  the  low  density  snow  near  the  dense  layer.  While  this  process  is  very 
slow,  it  can  produce  substantial  reductions  in  strength  of  the  low  density  snow  in  arctic  and 
antarctic  firn. 

Fabric  Gradient  Metamorphism: 

In  this  case  a  fine  grained  (high  energy)  layer  is  sandwiched  between  two  large 
grained  (low  energy  layers  of  the  same  density,  as  shown  in  Figure  12.  Figure  13  shows 
the  slow  migration  of  vapor  out  of  the  fine  grained  snow  into  the  courser  grained  layers. 
This  process  produces  slower  results  than  the  other  two  processes,  and  tends  to  be 
overpowered  by  the  other  two. 
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Figure  12.  Schematic  of  snow  cover  with  a  fine-grained  layer  sandwiched  between 
two  layers  of  course  grained  snow. 
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